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Abstract 

For  a  class  of  single-input  single-output  nonlinear  systems  with  unknown  con¬ 
stant  parameters,  we  present  a  direct  model-reference  adaptive  control  scheme, 
which  requires  only  output,  rather  than  full-state,  measurement.  The  nonlin¬ 
earities  are  not  required  to  satisfy  any  growth  conditions.  The  assumptions  on 
the  linear  part  of  the  nonlinear  system  are  the  same  as  in  the  standard  adaptive 
control  problem  for  linear  systems,  which  now  appears  as  a  special  case  of  the 
nonlinear  problem  solved  in  this  paper. 
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1  Introduction 


Until  a  few  years  ago,  adaptive  linear  [1,2]  and  geometric  nonlinear  [3,4]  methods  belonged 
to  two  separate  areas  of  control  theory.  They  were  helpful  in  the  design  of  controllers  for 
plants  containing  either  unknown  parameters  or  known  noniincaiities,  but  not  both.  In  the 
last  few  years  the  problem  of  adaptive  nonlinear  control  was  formulated  to  deal  with  the 
central  of  plants  containing  both  unknown  parameters  and  known  nonlinearities.  A  realistic 
plan  of  attack  to  this  challenging  new  problem  led  through  a  series  of  simpler  problems,  each 
formulated  under  certain  restrictive  assumptions.  The  two  most  common  assumptions  are 
those  of  linear  parametrization  [5-17]  and  full-state  feedback  [5-15]. 

The  purpose  of  this  paper  is  to  avoid  the  full-state  feedback  assumption  and  to  remove 
the  specific  restrictions  of  previous  output-feedback  results  [16,17]. 

In  the  linear  case,  the  adaptive  output-feedback  designs  follow  either  a  direct  model- 
reference  path  or  an  indirect  path  via  adaptive  observers.  Current  research  on  adaptive 
observers  for  nonlinear  systems  [18-20]  indicates  that  the  indirect  path  may  become  promis¬ 
ing  for  adaptive  nonlinear  control.  However,  the  major  stumbling  block  along  this  path 
continues  to  be  its  linear-like  proof  of  stability  which  imposes  restrictive  conic  conditions  on 
the  nonlinearities  [16,17],  Under  such  linear  growth  constraints  the  actual  nonlinear  problem 
is,  in  fact,  not  addressed. 

In  this  paper  we  formulate  and  solve  a  truly  nonlinear  output-feedback  problem  by  fol¬ 
lowing  the  direct  model-reference  path  of  Feuer  and  Morse  [21].  In  contrast  to  other  more 
popular  adaptive  linear  control  methods  [1,2],  the  method  of  Feuer  and  Morse  offers  a  possi¬ 
bility  to  prove  stability  without  any  growth  constraints.  In  a  companion  paper  [22]  we  have 
exploited  this  possibility  to  solve  a  full-state-feedback  adaptive  nonlinear  control  problem. 

In  this  paper  we  present  an  adaptive  output-feedback  result  without  nonlinearity  growth 
constraints. 

The  results  of  this  paper  apply  to  nonlinear  input-output  models  consisting  of  a  linear 
transfer  function  and  output-dependent  nonlinearities.  The  coefficients  of  the  transfer  func¬ 
tion  and  the  parameters  multiplying  the  nonlinearities  are  unknown.  For  th~  'inear  part,  the - 
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assumptions  of  minimum  phase  and  known  sign  of  the  high-frequency  gain  are  the  same  as 
in  the  adaptive  linear,  control  theory,  which  now  appears  as  a  special  case  of  the  nonlinear 
theory  presented  iij  t^his  paper. 

For  easier  understanding,  the  new  adaptive  scheme  is  first  designed  for  a  particular  system 
of  sufficient  complexity  to  be  illustrative  of  both  the  design  procedure  and  the  stability 
proper'  ?s  of  the  resulting  closed-loop  adaptive  system.  In  Section  2  we  design  the  adaptive 
scheme  for  this  system  and  then  prove  the  stability  and  tracking  properties  of  the  resulting 
adaptive  system  in  full  detail.  The  design  procedure  for  the  general  case  is  presented  in 
Section  3,  and  the  proof  of  stability  and  tracking  is  given  in  Section  4. 

Nonlinear  input-output  models  are  intimately  tied  to  state-space  equations  which  orig¬ 
inate  from  nonlinear  physical  laws  expressed  in  specific  state  coordinates.  In  Section  5  we 
give  a  state-space  form  of  the  class  of  nonlinear  plants  which  have  the  desired  input-output 
representation,  and  characterize  this  class  of  plants  via  a  set  of  geometric  conditions. 

2  Adaptive  Scheme  Design:  An  Example 

The  purpose  of  this  section  is  to  make  both  the  proposed  adaptive  scheme  and  the  main 
features  of  the  Feuer-Morse  method  more  easily  accessible  to  the  reader  with  the  usual 
background  in  control  theory  and  limited  familiarity  with  adaptive  linear  control. 

2.1.  Nonlinear  system  properties.  The  nonlinear  system  is  assumed  to  be  minimum- 
phase  [3,  Chap.  4]  and  its  nonlinearities  depend  only  on  the  output  variable.  This  implies  that 
the  nonlinar  system  is  linearizable  by  output  injection  [23].  The  input-output  description  of 
a  typical  nonlinear  system  of  this  kind  is  given  by 

D5y  =  {D2  +  2D  +  1  )u  +  6  D2p2(y)  +  Dp\{y)  +  po(y)  .  (2.1) 

where  u  and  y  are  the  scalar  control  and  output,  respectively,  D  —  and  0  is  an  unknown 
constant  parameter.  To  address  a  truly  nonlinear  problem,  we  choose  the  nonlinearities 
which  do  not  satisfy  linear  growth  constraints: 

Po(y)=y 3,  P\{y)  =  y2  F  2y3,  p2{y)  =  yey  +  2y2  +  .  (2.2) 
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It  is  important  to  notice  that  these  nonlinearities  are  not  in  the  span  of  u,  and,  hence,  the 
system  (2.1)  is  not  full-state  linearizable  by  static  output  feedback,  or  even  by  static  full- 
state  feedback,  as  shown  in  Sect.  5.  However,  it  is  input-output  linearizable  by  full-state 
feedback  [3,  Chap.  4]. 

The  above  structural  and  growth  properties  of  (2.1)  and  its  relative  degree  [3,  Chap.  4] 
show  that  (2.1 )  is  a  nonlinear  system  of  considerable  complexity.  However,  this  svstem  also 
satisfies  a  structural  constraint  under  which  the  results  of  this  paper  are  applicable:  the 
nonlinearities  do  not  enter  the  system  before  the  control  input  u. 


2,2.  Augmenting  the  CE  control.  As  in  most  adaptive  designs,  our  first  step  is  to 
find  a  dynamic  output-feedback  control  that  guarantees  the  specified  stability  and  tracking 
properties  when  the  parameter  0  is  known.  Most  adaptive  schemes  then,  replace  th  unknown 
0  with  its  estimate  0  and  implement  the  so  formed  "certainty-equivalence’'  control.  Such 
certainty-equivalence  designs  have  been  satisfactory  in  adaptive  linear  control,  but  have  failed 
to  produce  truly  nonlinear  results  because  of  their  inherent  linear  growth  constraints  [16,17]. 
To  avoid  this  difficulty  we  must  go  beyond  the  certainty-equivalence  approach.  Following 
Feuer  and  Morse  [21],  we  will  augment  the  certainty-equivalence  control  by  an  additive  term 
u  which  will  counteract  the  effects  of  rapidly  growing  nonlinearities.  It  will  also  provide  us 
with  additional  flexibility  in  the  proof  of  stability. 

The  certainty-equivalence  part  of  our  control  will  be  designed  to  match  a  reference  model 
of  the  same  relative  degree  as  that  of  the  nonlinear  plant  (2.1).  As  this  plant  is  input- 
output  linearizable  by  full-state  feedback,  we  will  choose  the  simplest  linear  reference  model 
of  relative  degree  three: 

( D  +  l)3yr  =  r  .  (2.3) 


The  first  step  in  matching  this  reference  model  is  to  filter  the  plant  equation  (2.1)  by  the 
strictly  proper  stable  filter  F  /  E2,  where  F  is  a  monic  polynomial  of  degree  2,  and  E2  is  a 
monic  Hurwitz  polynomial  of  degree  4.  This  results  in 


FA 

fJ 


FB 

E2 


u  4-  0 


'FD2  ,  ^  FD  ,  X  F  ,  X 

-rr-My)  +  — —pi(u)  +  —pu(y) 

Lj 2  £/  2 
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where  A  =  D5,  B  —  D2  +  2D  +  l  as  in  (2.1).  It  is  now  straightforward  to  verify  that  the 
desired  matching  is  achieved  by  the  control 


provided  that 


*  ^  ,  a  (  \  FB  -  E2  , 

u  =  —y  +  r  —  Ov\y) - - - u 

F  2  E2 


FD 2  FD  F 

v\y)  =  —^—Piiy)  +  -p-pi{y)  +  -rrpo(y)  • 

-C/2 


and  that  G,  a  polynomial  of  degree  4,  and  F  satisfy  the  polynomial  equation 

FD5  +  G  =  (D  +  1)3£2-  (2.7) 

Note  that  the  polynomial  FB  —  E2  in  (2.5)  is  of  degree  3,  since  FB  and  E2  are  both  monic 
polynomials  of  degree  4.  As  an  illustration,  the  choice  E2  =  {D  +  2)4  yields  the  following 
solution  of  (2.7): 

F  =  D2  +  11D  +  51  (2.S) 

G  =  129D4  +  192jD3  +  168£>2  +  SOD  +  16  .  (2.9) 

When  the  control  (2.5)  is  applied  to  the  system  (2.1)  and  the  initial  conditions  of  the 
filters  used  in  (2.3),  (2.5)  and  (2.6)  are  exactly  matched  with  those  of  the  system  (2.1).  then 
(2.5)  achieves  the  exact  tracking  y(t)  =  yr(t)  for  all  t  >  0.  However,  the  initial  conditions  of 
(2.1)  are  unknown  and  the  tracking  can  be  achieved  only  asymptotically,  that  is, 


y(t)  =  yr(t )  +  e(t)  -+  yT(t)  as  t  ->  oo , 


(2.10) 


where  e(t)  is  the  exponentially  decaying  tracking  error  caused  by  the  mismatch  of  the  initial 
conditions. 

When  the  parameter  0  is  unknown,  we  replace  it  in  (2.5)  by  its  estimate  0,  to  be  obtained 
from  a  parameter  update  law.  To  this  “certainty-equivalence”  part  of  our  control  we  add  a 
term  u  which  will  be  a  handy  tool  later.  So,  our  adaptive  control  will  be  of  the  toliowing 


form: 


G  -  .  ,  FB  -  E2 

u  =  y  +  r  —  0i/(y) - - - u  -f  u  . 

F  2  B2 


(2.11) 


When  applied  to  the  nonlinear  plant  (2.1),  this  control  yields  the  following  input-output 
description  of  the  resulting  feedback  system: 

y  =  (zTTIf  lr  +  (e~  uiy)  +  S1  +  ’  (2-12) 

where,  as  in  the  case  when  9  was  known,  e(t)  contains  all  the  exponentially  decaying  terms 
caused  by  the  mismatch  of  the  initial  conditions.  It  should  be  observed  that  with  an  exact 
estimate  9  —  9  the  linearization  of  (2.12)  is  achieved. 

Introducing  the  error  variables 

e  =  y-yr,  0  =  6-9,  (2.13) 

and  taking  the  difference  between  (2.12)  and  the  reference  model  (2.3),  we  obtain  the  tracking 
error  equation: 

e  =  (~pTT)3  +  t2]  +*(*)•  ('2-14) 

2,3.  Error  augmentation  and  swapping.  Following  the  standard  practice  in  adaptive 
control,  we  now  set  out  to  construct  an  error  equation  in  which  the  parameter  error  is  filtered 
only  by  a  strictly  positive  real  (SPR)  transfer  function.  As  a  first  step,  we  rewrite  (2.14)  in 
the  form 

e  =  coTTI  [VttfH  +  (077Ffi 

f (cTTF  M  "  WTT)  ['ioW'(!,)]  +  t(i) '  |215) 

The  first  term  in  (2.15)  is  in  the  desired  SPR  form,  while  the  second  term  is  due  to  the 
additional  control  term  u.  As  for  the  third  and  fourth  terms,  these  are  the  familiar  swapping 
terms ,  whose  presence  is  caused  by  the  time-varying  nature  of  0 :  if  0  were  constant,  these 
two  terms  would  cancel  out.  Let  us  therefore  define  the  augmented  error  e  as 


e  =  e  4-  7o  , 


(2.16) 


where  the  term  r/0  represents  all  the  undesirable  terms  in  (2.14): 

7/0  =  *  _(D+  1)3U  +  ( D  +  l)3  "  (D  4  1)  %D  + 


(2.17) 
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The  signal  multiplying  9  in  the  first  brackets  is  of  particular  importance  and  is  denoted  by 


h\  — 


1 


;v{y) 


C-US) 


{D  +  l)2 

Considering  v{y)  as  the  input  and  h\  as  the  output,  we  represent  (2.18)  in  the  state-space 


form 


where 


hi 

hn 


=  h  =  A0h  +  bu{y)  . 


,40  — 


0  1 

-1  -2 


b  = 


0 

1 


(2.1!); 


1 2.20  i 


It  can  now  be  verified  that  r;0  is  the  output  cf  the  third  order  system 


V  o  =  -Vo  +  Vi 


(2.21 


Vi 

V2 


=  V  =  -4or?  —  hu  —  hO 


i  2.22 ! 


The  variables  /ij  and  from  (2.19)  and  (2.22)  allow  us  to  express  the  tracking  erro>  as 

1 


e  = 


-  \0hi  -  7/xl  +  f\t) 


2.23) 


r  a 

The  analogous  expression  for  the  augmented  error  is 

1 


e  = 


£)-+T  \0h'\  +e(') 


(2.2-1) 


and  it  has  the  desired  SPR  form:  the  parameter  error  0  multiplied  by  the  "regressor"  h\  is 
the  input  into  the  SPR  filter  1/{D  +  1). 


2.4.  Update  law.  From  this  point  on,  the  route  prescribed  by  most  of  the  adaptive  linear 
control  literature  is  to  choose  a  normalized  gradient  update  law  and  to  set  u  =  0  (thus 
returning  to  a  pure  certainty  equivalence  control).  In  the  case  of  adaptive  linear  systems, 
boundedness  of  the  closed-loop  signals  can  then  be  established  using  the  Gronwall  lemma 
or  some  type  of  small-gain  argument.  Attempts  to  apply  this  type  of  stability  proof  to 
nonlinear  systems  have  so  far  been  successful  only  when  conic  constraints  are  imposed  on 
the  nonlinearities.  Without  such  linear  growth  constraints,  the  term  (^0  —  ())  n(;/)  can  cause 


i 


sonu-  signals  to  escape  to  infinity  in  finite  time  if  the  parameter  error  0  —  0  is  no!  rapid!) 
decreased.  1  lie  difficulty  with  normalizations  of  update  laws  is  that  they  don  t  allow  a 
rapid  enough  decrease  of  the  parameter  error  when  this  error  is  most  harmful.  A  simulation 
example  of  instability  of  a  full-state-feedback  scheme  with  normalization  [11;  is  given  in  our 
companion  paper  [22],  where  it  is  also  shown  that  an  unnormalized  update  law  preserves 
global  stability.  We  are,  therefore,  motivated  to  look  for  an  unnormalized  update  law 
I  lie  SPR  form  of  (2.24)  suggests  the  unnormalized  grade  nt  update  law 

0  =  hxr  .  .  2.25 , 


(liven  the  complexity  of  the  nonlinear  system  (2  1).  it  is  likely  that  such  a  simple  update 
law  will  shift  the  difficulties  in  the  adaptive  design  to  the  proof  of  stability.  Indeed,  this  is 
the  case.  A  simple  Lyapunov-like  function  involving  i2  and  (P  has  a  nonpositive  derivative, 
but  fails  to  prove  boundedness  of  y.  ft  clearly  shows,  however,  that  i/0  must  be  taken  into 
account.  Our  next  attempt  is  with  the  function 

V„  =  ^  (V  +  02  +  e2(r)drj  +  ~yx Pi]  ,  (  2.26  i 

where  P  is  the  positive  definite  solution  of  the  Lyapunov  equation 

PA0  +  AZP=  -2/  =►  P=  *  |  (2.27, 

L'sing  (2.22),  (2.24)  and  (2.25),  we  compute 

Vn  —  —  -  (e  —  e[t))2  —  -  je2  +  2rjX  Phh\'e  -f  2yT  Phiij  .  (2.28) 


and  try  to  render  it  nonpositive.  The  tool  we  have  prepared  for  this  task  is  it.  However,  it 
turns  out  to  be  impossible  to  counteract  the  effects  of  the  two-dimensional  vector  lihxr  by 
it  alone.  Hence,  we  need  an  additional  degree  of  freedom  involving  /i,.  This  prompts  us  to 
replace  tj  in  (2.26)  bv  the  new  vector 


Ci  -  'h 
C2  = 


pj.2fl) 
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fi' it! 


w*  ere  the  nonlinear  function  Cp/ql  is  at  our  disposal  With  the  new  wtriab|e>  2  i.h 
negative  function  to  he  used  in  our  proof  is 


f  =  -(<•*  +  +•  /  tz(T)<iT  |  +  - 

j  \  Jt  )  2 


C 1  fX 


lo  evaluate  l  we  need  (\  which  is  obtained  by  differentiating  (2.29)  and  using  (2.19;  '2.2'. 

Ci  =  Ci  -  M  ~  si^'i  )Cr 

sj  —  —  si  —  2c, 2  +  2C i ( ^ i  )C.  —  ii  —  h  >t)  i  2.9, 


+  si(M  C2  -  h,0  -  £.(/>.,  )Cj]  I  ~-/t.>Ci 

J  o'/?  j 


fnt  roducing  t  he  not  at  ion 


-Mo  -  Mi 


we  compute  the  flint'  derivative  of  l  as: 


f(M)2  -CTC  -  Uo  -  c'  /bet 


-ClP  c  -  fi.  M  ,  - 

-Cisi  +  a  —  s)W  +  MS  1  -  777— Ml 

f/n  1 


2.5.  Design  equation.  We  now  have  two  tools  to  make  V  nonpositive:  the  function  W./'i  . 
and  the  control  term  u.  With  these  tools  we  will  attempt  to  represent  the  quantity  enclosed 
in  braces  in  (2.33)  as  the  sum  of  two  squares.  It.  turns  out  that  this  is  possible  to  achieve  by 
decomposing  P  as  / 3  —  l\  -f  P,  such  that  the  following  design  equation  holds: 


P\wwl  /'’iC  r  Piii'iL^  P2C,  =  P 


—  2slCl  T  It  —  si  S2  +  sfCl 


I  he  substitution  of  (2.31)  into  ,2.33)  yields  the  desired  form  for  l  : 


V  --  -  Pt))2  -  C'C  -  -  C T P 1  irJ  -  ~  CTP:‘‘j  <  0.  (2.35  1 

Our  task  is  now  to  find  I\,  P,  £\\h\)  and  u  which  satisfy  the  design  equation  (2.31  1 .  For 
the  example  considered  here,  the  systematic  procedure  of  Sect.  3  gives  the'  following  solution 
tor  !’,  and  /h: 


Substitution  of  ( 2  if > )  into  (2.34)  results  in 


2  hlC 


1S1 


{(h[h-2  4-  h\  +  Ci  h\ )-  —  — /* i  j  Ci  +  [( h  i  h2  +  h\  4-  si  h\  )*]  C2 

6  Ci 

— 2Ci  +  Ci  —  Ci  —  C1C2  +  u 

whnh  direct.lv  yields  the  following  solutions  tor  ^  and  u: 

Cl  -  2 h* 

a  =  [s/if/io  4-  2/i,  -  4/if  4-  (hxh2  4-  h\  +  2/if  )2]  (h 

4  [2/1 1  +  {hi  h-2  +  h.'i  +  2/lj;  j  C'2  • 

I  'sing  (2  22)  and  the  notation 

y,  -  S h*h-i  4-  2/i  1  +  ( 1  4  2/i? )( h lh2  +  h'\  +  2/i,  )2 

y'2  =  2/i j  -(-  ( h  1  h2  4"  d ;  +  2/ijJ)‘ 

a  ran  be  defined  in  terms  of  available  signals  as 

U  =  T?i  +  r>2'/2  • 

To  summarize,  the  complete  closed-loop  adaptive  system  is 


Plant: 

DPj 

Control: 


=  [D2  +  2D  +  \)u  +  0  D2p2(y)  +  Dp^y)  +  p0{y] 


G  .  ,  x  FB  -  E2 

u  =  -  —y  4 -  r  -  Ou{y) - - - u  +  W/i  +  4V)2 

L2  t2 


Update  law: 


0  =  /i,  e  ,  e  =  e  4  r/0  =  y  -  yr  4  r/0 


l'  lit  c 


:  r.s; 


/)  =  ,4,,/t  4-  hid//) 

//  =  ,40r/  -  /)(y,r/i  4-  ^  /'  1 c 

f/o  =  ~  '/a  +'/n 


(2  37) 

12.34 ) 

(2.39) 

(2.  10) 

1 2  1 1  ; 

(2.42) 


(2.13 
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where  yr  is  the  output  of  the  reference  model  (2.3)  and  u(y)  is  defined  in  (2.6). 


2.6.  Stability  and  tracking.  The  stability  and  tracking  properties  of  (2.43)  are  now 
established  using  the  nonnegative  function  V  from  (2.30),  whose  derivative,  given  in  (2.36).  is 
nonpositive.  Because  of  the  piecewise  continuity  of  r(t)  and  the  smoothness  of  the  nonlinear 
functions  p0,  /h > P'2,  the  solution  of  (2.43)  has  a  maximum  interval  of  definition,  which  we 
denote  by  [0,(f).  We  will  now  show  that  t(  —  oc. 

From  (2.30)  and  (2.35)  we  conclude  that  e,  £,  and  0  are  bounded  on  [0,  ff).  Solving  (2.29) 
for  r)  and  using  (2.38)  we  obtain 


Pi  =  Ci 

P2  =  Ca  -  2A}Ci 


(2.44) 


Thus,  the  boundedness  of  £  implies  that  tjx  is  bounded,  »\..ich,  in  turn,  implies  that  r/0  is 
bounded:  r/0  =  -;/0  -f  j)l .  Since  e  and  y0  are  bounded,  e  is  bounded:  e  =  e  —  r/0.  By  the 
boundedness  of  yr.  this  implies  that  y  is  bounded:  y  =  e  +  yr  Hence,  v(y)  is  bounded,  which 
means  that  h  is  bounded:  h  —  A0h  +bu{y).  (2.44),  The  boundedness  of  h  and  (  implies  that 
y  and  u  are  bounded  (cf.  (2.44)  and  (2.39)). 

This  does  not  yet  prove  that  u  is  bounded.  From  (2.11),  to  prove  the  boundedness  of  u 


we  only  need  to  show  that 


FD  -  Eo 


u  is  bounded.  Since 


FB-E2  . 


F  B  E2  .  ,  f 

we  can  express  - — - u  in  the  torm 


e2 


is  of  relative  degree  1. 


£2 


FB  -  E2 

E, 


E2 


(D+  l)3(D  +  A,) 


£3 


(D+  1) 


-U  + 


1 


D  +  A,  (D+  l)3 


+  e(0 


(2.45) 


where  ,\  1  is  a  positive  constant  and  £3,  £.(  are  polynomials  of  degree  2.  Now  (2.45)  clearly 

D‘ 

shows  that  u  is  bounded  if  — rr — —  u  is  bounded  tor  1.  =  0,1,2.  Since  //  is  bounded  and 


1 


(D  +  l)3 

the  plant  is  minimum  phase  and  of  relative  degree  3,  it  follows  that  „ 

1  6  ( D  +  1  )3 

Differentiating  (2.14)  and  substituting  —0i/[y)  +  u  from  (2.11),  we  obtain  for  i  =  1.2 


■u  is  bomnh'd. 


10  _ 


ir 


{0+  1) 
i) 


-  [0n(y)  +  u  +({t) 


(D+\  ):1 


FB 


C 


7  jfMy)  +  —  "  +  —  / 


n  -  r 


+  r(/| 


(2.4(0 


1! 


Using  h, 


D 


1-1 


(P  +  1 


—  v(y),  i  —  1,2,  from  (2.19)  and  rearranging  terms  in  (2.46)  we  get 


p* 

FB  ■ 

— 

D' 

G 

(P  + 1)3 

.  e2  u. 

—  t 

(P+  l)3 

It2,j  - 1 

P 


P  +  1 


h,0  +  e(  t ) 


(2.47) 


From  (2.24)  we  have 


e  = 


e  = 


e  —  rjo  =  —e  4-  9hx  +  Vo  —  Vi  +  e(l) 

—  e  +  Ohi  +  9h\  +  Vo  —  Vi  +  c(0 
=  e  -  9hx  -  h\e  +  0/i2  ~Vo  +  Vi~  Vi  +  h\e  +  e(t) 


(2 -4S) 


>.49) 


Since  e,  /i,  7/,  0  are  bounded,  (2.4S)  and  (2.49)  imply  that  e,  e  are  bounded.  Hence,  by  (2.47). 


D‘ 


rFB  1 
£2  U 


(D+  l)3 
follows  from  the  boundedness  of 


is  bounded  for  i  =  1,2.  The  boundedness  of 

1 


D‘ 


D' 


(D  +  l)3 
Dl 


(P  +  1)3 

u  and  the  recursive  expression 


u  for  i  =  1,2  then 


■u  = 


FB  D(FB  -  E2)  P-1 

u - 


:U  +  t{t) 


(2.50) 


(P+1)3"  (P+1)3P2“  E2  (P  +  1)3 

Next,  ve  prove  that  the  state  of  the  plant  is  bounded.  From  (2.12)  it  follows  that  D'y, 
0  <  i  <  3  are  bounded.  Combining  this  with  the  fact  that  the  plant  is  minimum  phase,  we 
conclude  that  the  state  of  any  minimal  realization  of  (2.1)  is  bounded  on  [0,  If). 

Thus,  we  have  shown  that  the  state  of  the  closed- loop  adaptive  system  (2.43)  is  bounded 
on  its  maximum  interval  of  existence  [0,  if).  Therefore,  If  =  oo. 

Finally,  we  prove  convergence  of  the  tracking  error  to  zero.  From  If  —  oo  and  (2.35)  we 
conclude  that  V  is  bounded  and  integrable  on  [0,oc).  Furthermore,  the  boundedness  of  e 
(cf.  (2.24)),  C  (cf.  (2.31)),  and  h  (cf.  (2.19)),  implies  that  V  is  bounded.  Hence.  \  — ►  0  as 
I  — +  oc,  which  implies  (cf.  (2.35))  that  e  — >  0,  (  — ►  0  (since  h  is  bounded).  This,  in  turn, 
implies  that  Vo  — *  0  as  /  — ♦  oc  by  (2.21).  We  conclude  that 


>J  ~  Ur  =  ('.  =  e  - Vo 


0  as  t  —*  oc  . 


[2.5 1 ' 


□ 


3  The  Systematic  Design  Procedure 


Even  though  the  expressions  in  the  general  case  become  more  complicated  than  in  the 
preceding  section,  the  main  steps  of  the  design  procedure  remain  essentially  the  same. 

3.1.  Nonlinear  system  properties.  We  consider  the  class  of  n-dimensional  nonlinear 
systems  which  have  an  input-output  description  expressed  globally  by  the  n-th  order  scalar 
differential  equation 

A(D)y  =  B{D)  [<7(j/)u]  +  ]T  D'  [pl0(z/)  +  pj  (y)01\  .  (3.1) 

i=0 

where 

•  the  coefficients  a0,  ■  .  ,an-i  of  the  denominator  polynomial  A{D)  =  Dn  +  an_lDn~l  + 

+  «o  arc  unknown, 

•  the  coefficients  60, . .  . ,  bm(m  <  n  —  1)  of  the  numerator  polynomial  B(D)  —  b,nDm  + 

■  +  b0  are  unknown,  but  B(D)  is  known  to  be  Hurwitz,  and  the  sign  of  bm  is  known. 

•  0 ,  is  an  ^-dimensional  vector  of  unknown  parameters, 

•  <7(r/ ),  Pij(y ),  0  <  i  <  m,  0  <  j  <  t,  are  smooth  nonlinearities  with  q(y)  ^  0  Vy  £  IR. 

Pij( 0)  =  0,  0  <  i  <  m,  0  <  j  <  i. 

Systems  in  this  class  are  linearizable  by  output  injection,  and  input-output  linearizable  by 
full-state  feedback,  but  not  necessarily  full-state  linearizable,  even  by  full-state  feedback,  as 
will  be  shown  in  Sect.  5. 

3.2.  Augmenting  the  CE  control.  The  design  objective  of  the  certainty-equivalence 
part  of  our  control  is  to  match  a  reference  model  of  the  same  relative  degree  as  that  of 
the  nonlinear  plant  (3.1).  As  this  plant  is  input-output  linearizable  by  output  feedback,  we 
choose  the  linear  reference  model: 


B\{D)E2{  D)yr  -  Hr,  E2(D)  =  E2l(D)E,2(  [)) , 


where  E\{D),  E2(D),  E2X{D ),  E22{D)  are  monic  Hurwitz  polynomials  of  degree  n  —  m.  n—  1. 
n  -  rn  —  1,  and  m,  respectively,  and  R{D)  is  a  polynomial  of  degree  h  <  n  —  1.  Filtering  (3.1 ) 
by  the  strictly  proper  stable  filter  F/  E2,  where  F  is  a  monic  polynomial  of  degree  n—m—  1. 
we  obtain 

FA  FB  r  ,  ™  FD‘  r  T  , 

-£-*/  =  —  [<i{y)u]  +  Y  ~~E'  [p.-o(y)  +  pj{y)8 i]  •  (3.3) 

2  2  1=0  2 

It  is  now  straightforward  to  verify  that  in  the  case  when  the  coefficients  of  A(D)  and  B(D ) 
and  the  parameters  9X  are  known,  the  desired  matching  is  achieved  by  the  control 

1  (  G  M  ,  J  Lo  f  ,  .  .,1 

provided  that  G'.  a  polynomial  of  degree  n  —  1,  F,  iV/,  T:  and  L0  satisfy  the  polynomial 


equations 


FA  +  bmG  =  EXE2  (3.5) 

M  =  ^-R  (3.6) 

Tj  =  ~0ijF,  0  <  j  <  ^ ,  ^io  =  1  (3.7) 

Go  =  ^-FB-E2.  (3.8) 

Note  that  Lo  is  a  polynomial  of  degree  n  — 2,  since  both  7 —FB  and  E2  are  monic  polynomials 
of  degree  n  —  1.  When  the  control  (3.4)  is  applied  to  the  system  (3.1).  asymptotic  tracking 


is  achieved: 


1/(0  =  yr (0  +  e(0  — ►  yT(t)  as  t  -*  00  , 


where  f.(t)  is  the  exponentially  decaying  tracking  error  caused  by  the  mismatch  of  the  initia 
conditions 


We  now  rewrite  the  control  (3.4)  as 


= - T-T  u’,r)0 

<i(y)  1 


(3.10) 


where  the  tig dimensional  vectors  0  and  tj>  (with  rig  =  2n  +  (f.  +  l)(n  —  rn))  are  defined  as 
9  [<70,  ■  .  fjn  —  1  1  n>() '  Ll(t ,  ■  ■  •  ,  — m  —  1  >  G O,  ■  ■  .  G  n  —  rn  —  1  ' 


....  Go,  •  -  ,  tafn-m-  1,  G)0-  •  •  ■  ,  G)n-2j 
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o  (y,  u\  r ) 


Dn_,  R  Dn_m_!  D'  ,  x 

~Wy'%r' 


p  ^  p  PM '  p  Ps  p 

&  21  i=0^27  &2\  x=QL22 


T.~Er~P'l{,j) 


D 


n  —  m  —  1 


.  .  £>„_2 


Ht —P,e{y), 


£21  ,=0^22 


[<?(»/)“*] 


(3.12) 


with  0^  defined  as  the  (A:  +  l)-dirnensional  row  operator  Dk  =  [1 ,  D, .  . . ,  Dfc].  The  form 
(3.10)  is  particularly  useful  in  the  case  where  the  coefficients  of  A(D)  and  B(D)  and  the 
components  of  the  vector  0y  are  unknown.  Since  in  that  case  the  parameter  vector  o  defined 
in  (3.11)  and  used  in  (3.10)  cannot  be  computed,  it  is  replaced  by  an  estimate  0.  The  so 
formed  “certainty-equivalence”  control  is  then  augmented  by  an  additive  term  u  which  is  yet 
to  be  designed.  Hence,  the  adaptive  control  will  be  of  the  form: 


1 


u  = 


q(y) 


-(pT(y,u,r)0  +  u 


Filtering  the  system  equation  (3.1)  by  the  strictly  proper  stable  filter 


F 


EyE2 


(3.5)-(3.S)  and  (3  11)-(3.12),  we  obtain: 


FA 

EyE 


■y  = 


BF 

EyE2 
/ 
by 

Ey 

6, 


(?(»)“]  +  E  ppr  \pio(v)  +  +  £(<) 


y  =  T-  { jh2  W»W  +  §;y  +  E 


T  m  D’ 

_  .jr-ILy-Mv) 

j=0  1*21  j_o  *22 


+  e{t) 


(3.13) 
and  using 

(3.14) 

(3.15) 

(3.16) 


y  =  |<?(2/)w  +  <t>T{y,  u,  r)0  +  ‘—T^  +  e(t), 

where,  as  in  the  case  when  0  is  known,  e{t)  denotes  a  linear  combination  of  exponentially 
decaying  terms  caused  by  the  mismatch  of  the  initial  conditions.  Substitution  of  (3.13)  into 
(3.16)  yields  the  following  description  of  the  resulting  feedback  system: 


y 


R 


Ex  \bmE2 

Introducing  the  error  variables 


+  0T(y,u,r)  [d  -  0)  4-  iij  +  e(t) . 


e  =  y  -  J/r ,  0  =  0-0 


(3.17) 


,3.181 


and  taking  the  difference  between  (3.17)  and  the  reference  model  (3.2).  we  obtain  the  tracking 
error  equation: 


Ey 


<pl(y,u,r)0  +  u  +  f(t) 


3.1!)) 
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In  the  special  case  of  relative  degree  one  (n  -  m  —  1),  the  design  is  extremely  simple. 
Since  the  transfer  function  |6m|  f  Ex  is  SPR,  the  parameter  update  law 

6  =  sgn(6m)r</»(j/,u,r)e,  (3.20) 

where  T  =  PT  >  0  is  the  adaptive  gain ,  guarantees  boundedness  of  all  the  closed-loop  signals 
and  convergence  of  the  tracking  error  e  to  zero  [1,  Chap.  5],  and  the  control  augmentation  is 
not  needed:  0  =  0. 


3.3.  Error  augmentation  and  swapping.  For  relative  degree  higher  than  one  ( n  —  m  >  1 ). 
the  design  becomes  considerably  more  complicated,  since  |6m|/£’1  is  no  longer  SPR.  We  first 
rewrite  (3.19)  in  the  form 


e  = 


1 


1 


bm9T—<t>{y,u,r) 
& o 


D  -f-  Aq 


b-m  - 

+  —  U 


D  +  Aq  [Eq 


Ex 

+  e(<) , 


(3.211 


where 


Eo(D){D  +  A0)  =  EX{D) .  (3.22) 

In  contrast  to  the  example  of  Sect.  2,  where  the  high-frequency  gain  bm  was  known,  here 
it  is  unknown.  Therefore,  using  an  estimate  bm  and  denoting  6m  =  bm  —bm,  we  rewrite  (3.21 ) 


in  the  form 


1 


e  = 


D  - f  A0 
1 

+  ■ 


b. 


i9Tjp-<t>  +  ~bm  ( -pru  -  0T-jp-<t>  + 
tjQ  \  Eq  Eq  Eq 


0To 


,  ,  1  rT  1  1 

bm  !  -7 r«  -  e  T -<t>  +  TT 

Eq  Eq  Eq 


drd> 


+  e(0  • 


(3.23) 


D  +  A0 

Since  the  first  summand  in  (3.23)  is  in  the  desired  SPR  form,  we  define  the  augmented  error 


e  =  e  +  g  o  , 


(3.2-1) 


where  the  term  r/0  represents  all  the  undesirable  terms  in  (3.23): 


1 


m  =  - 


D  +  A0 


,  /  1  -  /IT  1  1 

bm  (  ~p~u  ~  +  — 


0V6 


(3.2a) 
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The  vector  multiplying  6  in  the  first  term  is  denoted  by 

1 

0  =  —<p{y,u,r)  . 

e-0 


(3.26) 


Considering  <p(y,u,r)  as  the  input  and  0  as  the  output,  we  represent  (3.26)  in  the  matrix 
state-space  form: 


H  =  AqH  +  b<pl{y,u,r) 


(3.27) 


0T  =  CT//  , 

where  (c,  ,40,  6)  is  a  minimal  realization  of  1  /Co : 

c?{sl  —  Ao)~lb  = 


EoU 


(3.28) 


Now  7o  1S  the  output  of  the  ( n  —  m)-dimensional  system 

t)o  =  —  +  bmC^T) 

i)  =  AoV  —  bu  —  HO 


(3.29) 

(3.30) 


The  variables  0  and  cTr]  from  (3.27)  and  (3.30)  allow  us  to  express  the  tracking  error  as  (cf. 
(3.21)) 


1 


e  = 


6m0T0  -  bmCTT)  +  e(t) . 


s  +  A0 

The  analogous  expression  for  the  augmented  error  is 

1 


(3.31 : 


s  +  Ar 


bm0lxl>  -  bmc  Tj\  4-  f{i) , 


H.32' 


and  it  has  the  desired  SPR  form:  the  parameter  errors  0  and  bm  are  filtered  only  by  the  SPR 
filter  l/(s  +  A0). 

3.4.  Update  law.  As  in  Sect.  2,  we  choose  the  unnormalized  gradient  update  laws  suggested 
by  the  SPR  form  of  (3.32): 


0  =  sgn(6m)r  0e 

bm  =  ~ ~fCTr]e , 


(3.33) 

(3.31) 


where  I'  =  fT  >  0  and  7  >  0  are  the  adaptive  gains. 


From  Sect  2  vve  knew  that  in  the  proof  of  stability  there  will  be  a  need  to  balance  t he 
interaction  between  t]  and  H .  Therefore,  we  introduce  the  new  variables  (  : 


(=S~lT),  S  =  C\l=.Cn, 


(3.35 ) 


where  n  =  n  —  m  —  1. 


C,  = 


c,Ajc,...,  (aJ)‘  c 


lT 


,  i  =  1 , . . . ,  h 


(3.36) 


— '  —  ^nxn  T 


0 

i2E2.fi 


(3.37) 


in  —  1  En-l.n 

Ej.i  =  ’  0  <  j  5:  *  1  (3-3S) 

and  /txt  is  the  1  x  1  identity  matrix.  The  components  of  the  z-dimensional  row  vectors 
are  nonlinear  functions  of  the  elements  of  H  which  represent  the  aforementioned  additional 
degrees  of  design  freedom.  In  order  to  show  that  the  matrix  S  defined  in  (3.35)  is  invertible, 
we  note  that,  because  of  the  structure  of  the  matrices  defined  by  (3.38),  the  matrix  E 
is  lower  triangular  with  ones  on  its  diagonal.  From  this  and  the  aforementioned  functional 
dependence  of  i,  on  H ,  it  follows  that  Z-1  always  exists  and  that  the  elements  of  both  E 
and  E"1  are  polynomial  functions  of  the  elements  of  H  Furthermore,  Cj"1  exists  because 
(cA,  /f0)  is  assumed  to  be  an  observable  pair. 

The  nonnegative  function  to  be  used  in  our  stability  proof  is 


v  =  \  I  K  I  0Tr-'i  +  Ul  +  j-  jT  e2(r )*■)  +  P(  . 


(3.39) 


The  form  of  (3.39)  is  the  same  as  that  of  (2.30),  where  P  is  the  positive  definite  solution  of 
the  Lyapunov  equation 


P  ,40  +  /10  P  —  —Q 0  • 


(3,10) 


To  evaluate  V  we  need  (,  which  is  obtained  by  differentiating  (3.35)  and  using  (3.30),  (3.33) 
and  (3.36)— (3.38): 
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C  =  s  1  ( AqT 1  -  bu  -  sgn (bm)HTipe)  +  (S  ')  t) 

=  S-lAoSC  +  S~l(bu  +  sgn  (bm)H  Pipe)  -  S~lSC 

=  A0C  —  sgn(bm)S~1  HTipe S~l(A0S  —  SA0  -  S)C  —  S~1bu  .  (3.41) 


Introducing  the  notation 


we  rewrite  (3.41)  as 


w  =  -sgn (bm)S  lHF ip 

W  =  S_1  (SA0  -  +  S) 

(  =  AqQ  +  we  —  W(  —  S~l  bu , 


(3.42) 

(3.43) 

(3.44) 


and  compute  the  time  derivative  of  V  as 

2 


C  =  - j  (A*  -  7^)  -  TT^oC  -  V  { ^  -  ?p  («*  -  -  5- to) 


(3.45) 


3.5.  Design  equation.  The  tools  we  have  at  our  disposal  to  make  V  nonpositive  are  the 
functions  £,(//)  and  the  control  term  u.  With  these  tools  we  will  attempt  to  represent  the 
quantity  enclosed  in  braces  in  (3.45)  as  the  sum  of  n  squares.  It  turns  out  that  this  is  possible 
to  achieve  by  decomposing  P  as  P  —  5K,n=i  Px  such  that  the  following  design  equation  holds: 


PtwwT PC  =  P{W C  +  S~lbu) . 


( 3.4G) 


t=i 


The  substitution  of  (3.46)  into  (3.45)  yields  the  desired  form  for  V: 

2 


1 


1 


2A0  f  e 


1=1 


<  o. 


(3-47) 


Our  task  now  is  to  find  P,,  ( II )  and  u  which  satisfy  the  design  equation  (3.46).  Following 

the  development  in  [21],  we  define 


P, 


CjMiCx . 


(3. IS) 
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where 


A/>  =  (C.P-'C^1 

Mi  =  (C.P-1^)"1  ^/lX, -'pCiP~lCjMiEj)j,  *  =  2, . . . ,  n  . 
Tn  [21.  Lemma  1]  it  is  proved  that  (3.48)-(3.49)  result  in 

Y2p,  =  p 

«=i 

CjP-'CjMiCi  =  0 ,  1  <j<i<n. 


This  proof  is  now  given  for  completeness.  From  (3.49)  we  have 

/.x.  =  CxP~xCj  Mi  +  CtP_iCj  MjEj'i 

j= i 

=  CiP~li2CjMk  Ek<i. 

t=i 


(3.49) 


(3.51 ) 


(3.52) 


Cfc  =  £fc,,C, ,  L-  <  j  ,  (3.53) 

we  obtain 

c,  =  C.P-1  £  CjT/fcCfc .  (3.54) 

it=i 

Evaluating  (3.54)  at  z  =  n  and  using  the  nonsingularity  of  Cn  we  obtain  (3.50).  Furthermore, 
premultiplying  (3.54)  by  where  j  <  i,  and  using  (3.53)  again,  we  obtain 

Cj  =  C2P-X  J2  CjMkCk  ,  i>j.  (3.55) 

*=i 

But  from  (3.54)  wc  have 

Cj  =  CjP~l  ^2  CjMkCk  , 

k=  1 

which,  combined  with  (3.55),  results  in  (3.51). 

Having  established  (3.50)  and  (3.51),  we  now  set  out  to  find  £,(//)  and  it  which,  along 
with  P,  defined  bv  (3.48)-(3.49),  satisfy  the  design  equation  (3.46).  Substituting  (3.48)  into 
(3.46)  we  obtain 

WC  +  S~lbu  =  P~x  £  Cj  Mowin'1  CjMjC, .  (3.56) 

i=i 
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Using  (3.43)  we  rewrite  the  design  equation  (3.56)  as 


(SA0  -  A0S  +  S)(  +  bu=  SP~l  £  Cj  M:n.ivtvTC?M?CA  . 


(3.57) 


i  =  1 


Premultiplying  both  sides  of  (3.57)  by  ct/4q  1  for  i  =  1, . . . ,  n,  we  obtain 


cT Al0~l{A0S(  -  bu)  =  cMq  1 


5  40-  P"1  CJMjCjWu^CJMjCj  +  S 


j=i 


+cT  Aq  lbu  ,  1  <  i  <  n  . 


(3  58) 


From  (3.35),  (3.37)  and  (3.53)  we  have 


CnS  =  ZCn  =  C*  + 


0 

Sn- 1(-  h-i 


which  gives 


(3.5!)) 


cTS  = 


ctA{,-,S 


=  cT A'q  1  +  £t_iC,_ i,  1  S  !  S  n 


(3.60) 

(3.61) 


Furthermore,  from  (3.35)  and  (3.37)  we  have 


cJA‘0-lS  =  cTA'0-lC^±Cn  = 


cTA\rlC;1 


0 

fiCi 


L  (a-lCfi-l 
£i  —  1  C i  —  1  i  1  —  *  —  H  ! 


where  we  have  used  the  definition  of  Cft  to  obtain  the  last  equality.  Finally,  from  the 
definition  (3.28)  of  the  tr.ple  (c,  A0,  b)  we  have: 


ct/1q  '6  =  0  ,  [  <  i  <  h  — 


zrA*~lb 


,3.6  1) 
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Substituting  (3.60 )— ( 3.64 )  into  (.3.58)  results  in 


CT,t‘  +  ZiC.  =  +  C1  .-1q  + 

n 

X  p-^C] MjC}wwtcJ UjCj .  1  <  i  <  n  —  l 


f  3.1).') 


j=i 


—  \r^  4^S’  -  C  C  —  rL  4*  _  C  (  '  A  i  I  r  1  4 

—  I  ^  s  a  —  1  -  n  —  1  £  -  sn-P-  ri-  !•  H)  •  •  * 


T  i  u  -  1 

U 


C  (  ’ 
<,n-  1  *  -  1 


X  p-lY cj .\1,C,W a;T c J M?C. 

(  *  J  •»  J  J  j  J 


<,  • 


3.6(1 


J=‘  J 

At  this  point,  we  have  almost  achieved  our  goal  of  findm  S  s.  and  u  which  satisfy  the  design 
equation  and  thus  render  V  nonpositive.  Still,  ( 3.65 )— ( 3.66 )  are  in  a  rather  complicated  hum 
and,  moreover,  they  involve  the  time  derivatives  of  the  functions  Therefore,  we  now  -•<-! 


out  to  simplify  (3.65)-(3.66)  and  to  express  ,  <as  explicit  functions  ot  available  signals. 

Motivated  by  the  appearance  of  the  terms  C},r  in  (3.65) -( 3.66 ),  we  introduce  t lie  ;■ 
dimensional  column  vectors  uq.  .  .  ivh  which  are  defined  as 


il\  —  C,  w  ,  1  <  i  <  n . 


i  3 .  t » i 


Combining  (3.67)  with  (3.53)  we  see  that  these  vectors  satisfy  the  recursive  expression.' 

a1*  =  C\w  TbSi 

w,  =  E,.l+lwl+[  ,  1  <  i  <  n  —  1  .  « 3.6f *  * 


Using  (3.38)  we  can  rewrite  (3.69)  as 

wj  =  u’i,i]  ,  1  <  i  <  n  —  1  .  (  3. 1 1)  i 

We  now  set  out  to  obtain  explicit  expressions  for  aq,..  ,  tr’a  in  terms  of  f,.  Sulxt  it  ut  uni 
(3.42)  into  (3.68)  results  in 

Hu;ft  =  -sgn(/jm)Cii//rC  <3  71  ■ 

We  then  use  (3.27).  (3  37)  and  (3.70)  to  rewrite  (3.71)  as 


/ 

0 

\ 

[  r 1  1 

Si  Ul 

ICr i  1 

=  -sgn  (bm) 

cTTo 

I  ii  X  n  3" 

n  —  1 

.71 

V 

c  /n 

.  Sa-lljn-l,n  . 

1 

.  r  1  .‘h"~  *  . 

IK  (d.d.S).  (3.72)  is  e»pn valent  to 


A  n  -  I  I  x  I  a  -  1 )  A 


u  \ 

Si  t-  l.fi-  1 

tc  n  —  1 

Sn  —  l  n  -  2.n  —  1  .  / 

^  h.ri  ”t~  6"i  —  l^n  —  1 


-sgn(  6m  )C\i_  (  //  r  II  r 


-sgn(6m)cI.*tr,//r//If 


Starting  from  (di’di  (-5  71).  one  can  repeat  the  above  procedure  to  show  that  i.  '3.  >  2 1  i> 


equivalent  to 


>‘~2J  +  si  ,('i 


=  -s%n(bn 


II VH 


L  +  ,  J  l  ^r1  J 

Hence  the  explicit  expressions  for  the  vectors  u  ,.  .  .  wH  are 

ii’l  -  -sgnl  h,„  )  ^cr  II I  II 1  c) 

»\r  -  Ulv  -sgn|/>m)(cr,i;r1//r//rc)  -6-.u-.-ij  •  '  = 

( 'ombinmg  id  ol )  and  (d.(>7)  we  obtain 

c\  '/  \ijC}u'u-Tc'!  mJ  Cj  =  r.p-^clu^w^cjyijc, 

j - 1  r=i 

=  (\  (\ 
J  =  1 

=  r,.v,c,. 


])  .  \  ■ 


\vll<Tr 


v.  =  r)-iZ<'jMJt>yrJuJi-:J,l.  i  =  i...  . 


I  Ini',  the  last  term  in  (d  bo)  and  (d  (>b)  can  lie  rewritten  as 


•  c1  = 


=  f[0.  0  1]  + 

j  =  t 

=  f|o.  in;  +■  6-i )(\\c, 

=  (rl  a;-'  +6-H’.-i).v/’. 
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Introducing  the  i  x  (i  +  1)  matrix  /i\  +  1  =  [0,  /,Xlj  and  Substituting  (3.80)  into  (3.65)(3.66) 
we  obtain 


Slf  l  —  S!-lf  t-1  +  Sl-lC 

\_i-4o  -  (cT^_1  +  ^-iC,_1).V,C’1 

C  J?  1  C 

1 S  i  -  l  ^  t  —  1 ,  t  i  Si  — 

,/?,  -(c'S4‘-1+^_1C,_1)rV1]C1.  i=  1,.. 

.  .  n  —  1  (3.81 ) 

\L  =  J 

{‘JX$  -  c1  .  l," 

|sn—  t  E'fi—  1  ,fi  —  \  Rfi 

-  (mg-1 

+  su-iCn-i)  -S]  cft|  c  ■ 

( 3.83 ) 

Hus  tonn  makes  it  apparent  that  tiie  design  equation  (3.56)  is  satisfied  bv  the  recursive 
expressions 

Sf.  =  -cT.V,  (3.83) 

si  =  +  6_i/?,  -  (cT^-'  +^,_1C)_i)  .V,.  i  =  2, .  . . ,  n  .  1 3.M ) 

«  =  [cr.40n  ~  (cT.4*  +  ^Cn)  S'1 1  r, .  ,  3.85) 

l’o  finally  solve  the  design  equation,  we  need  to  express  in  (3.84)  as  an  explicit  function 
of  available  signals  We  frst  show  via  an  induction  argument  that  (3.76).  (3.77).  (3.7'di. 
i  4  S3)  and  (3.84)  imply  that  the  elements  of  tn,,  .V,  and  are  polynomial  functions  of  the 
elements  of  C,  II : 

•  For  i=l.  this  fact  is  obvious  from  the  definitions: 

"d  =  —  sgn(6m  )C'i 

•V,  =  P-'Cj  Myw.wjMj 

SC1  =  -cT  V i 


for  i  —  k  <  n  —  1,  suppose  the  elements  of  tet,  .Vfc  and  are  polynomial  functions  of 
the  elements  if  f •*.-//.  I  hen.  the  derivative  of  ^  can  be  expressed  as 


s  k 


y'/j  t 


(3.86) 


where  /q.  is  the  j-th  column  of  CkIl  and  o(,k/<)hk,j  the  k  x  j  matrix  of  partial 
derivatives  of  fc- K- *.-//)  with  respect,  to  the  ehunents  of  /i*. }  But  from  (3.37).  (3.63). 
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(3.64),  and  the  definition  of  Rk+ , ,  we  obtain 


C\H  =  CkA0II  =  R-k+ 1 C h  1 1 H  ,  1  <  k  <  n  -  1 


(3.87; 


Combining  (3.86)  and  (3.87)  we  can  express  £k  as 


"8 


fc  =  Etf+i^+. 

j=i 


_^L 


(3.88) 


•  For  i  =  k  +  1,  we  have 


wI+i  =  [w?,  -sgn(6m)[0  ...  0  l]CA:+1//r(C’1//)T  -  £kwk] 

=  NkEk.k+l  +  F  1  Cj+lMk+lU>k+lwJ+iMk+l 

s*+i  =  E  ^4-1  .^Li  -£<a-+i  +  £,tRk+i  -  (c^/fp  +  ^fcCfc)Arfc+i  . 
j  =  l  vhk,] 

Hence,  the  elements  of  iwjt+i,  iV*+I  and  6t+i  are  polynomial  functions  of  the  elements 
of  Ck+ i H 

Thus,  the  term  t  in  (3.84)  can  be  calculated  explicitly  from  (3.88). 

The  design  procedure  is  now  complete  The  expressions  for  £,  and  u,  which  guarantee 
that  the  nonnegative  function  C  in  (3.39)  has  the  nonpositive  derivative  (3.47).  are 

,  h 

(3.89) 

The  designed  closed- loop  adaptive  system  is: 


c 

si 


c 

Si 


Tp 


u 

T 


-cT.V, 

V'  lT  oT  ^6-1 

E  Kjr  .  or — 

j= i  ani-\,j 

■plv 


<'TX  -  (cTd o‘  +  4nCn)  (C-'Ea) 


+  6-1 -  (cT.do-’  +  6-iC-i)  -V.  •  i  =  2 


-1 


‘) 


) 


Plant: 


i  =  Arx  +  brq(y)u  +  Y^-n-x  Pxo(y)  +  pj{y)0\ 


Control: 


u 


Update  law: 


-7—  f-0T(y,u,r)^  +  Yv 
<7  id)  1 


8  =  sgn(6m)r;/’e 
bm  =  -7  cT7?e 


e  =  e  +  rjo  =  y  -  yr  +  no 

Filters: 

H  —  A0H  +  boT  (t/,u,r) 
t/>T  =  cTf/ 
fio  =  —^oPo  +  bmcF  r/ 
r /  =  A)?;  -  VTr;  -  //sgn(6m)r0e , 


where’  ,//r  is  the’  output  of  the  reference  mode!  (3.2),  o(y.u,r)  is  defined  in  (3.12 
(cv,  /U-,  Ik  )  is  a  minimal  state  representation  of  the  plant  equation  (3.1): 


0 


-«n-l 

'  1  ' 

I 

0 

1  cr  = 

0 

bm 

-«o 

0..0. 

.  0  . 

L  J 

with  the  ( n  -  t)-th  coordinate  vector  in  Hln . 

The  stability  and  tracking  properties  of  (3,90)  are  established  in  the  ne’xt  section 


(3.90) 


).  and 


(3.91) 


26 


4  Stability  and  Tracking 

We  are  now  ready  to  state  and  prove  our  main  result: 

Theorem  4.1.  For  any  uniformly  bounded  and  piecewise  continuous  reference  input  r.  all 
the  signals  in  the  closed-loop  adaptive  system  (3.90)  are  well-defined  and  uniformly  bounded 
on  [O.oc),  and.  in  addition. 


lim  e(t)  =  0,  lirn  Tj(t)  — *  0,  lim  T}0(t)  =  0  .  (4.1) 

t — ►oo  t—*  oo  t — »oo 

Proof.  Due  to  the  piecewise  continuity  of  r(t)  and  the  smoothness  of  the  nonlinear  functions 
appearing  in  the  definitions  of  various  terms  in  the  closed-loop  system  (3.90),  the  solution 
of  (3.90)  has  a  maximum  interval  of  existence  [0,  tf).  On  this  interval,  the  time  derivative  ot 
the  nonnegative  function  V  defined  in  (3.39) 


v  =  i(*>+iMrr-i  +  i  + 


J-  f  t2iT)dT]  +  ~(TP(  , 

Aq  J  t  In 


computed  along  the  solutions  of  (3.90),  is  given  by  (3.47): 


We  conclude  that  V .  e,  6,  ~bm  and  (  are  bounded  on  [0.  If)  by  constants  depending  only  on  the 
initial  conditions  of  (3.90).  This  implies  that  0 ,  bm  are  bounded  on  [0,tf).  The  boundedness 
of  f  together  with  (3  60)  implies  that  cTr;  is  bounded;  from  the  definition  of  T]0  in  (3.90)  and 
the  boundedness  of  bm  and  cTr?  it  follows  that  1S  bounded.  But  since  e  =  e  —  rj0.  and  c.  q0 
are  bounded  we  have  that  e  is  bounded.  Now  from  the  boundedness  of  r  we  have  that  yr 
is  bounded,  and,  hence,  y  is  bounded,  since  ij  =  e  +  yr.  The  boundedness  of  y  implies  that 
all  the  nonlinearities  appearing  are  bounded,  and,  furthermore,  that  q(y)  is  bounded  away 
from  zero.  Filtering  the  system  equation  (3.1)  with  the  strictly  proper  stable  filter  l  /  BF\ 
and  rearranging  terms,  we  obtain 


1  rn  JJt 

Y] ■  +  1Tex  +  (y)°l 


which,  by  the  boundedness  of  y,  implies  that  ~q{y)u  is  bounded.  The  boundedness  of  H  is 

E\ 

now  established  by  proving  that  the  row  vectors  cT A'0H,  0  <  i  <  h  —  1  are  bounded.  This, 
in  turn,  is  proved  as  follows:  First,  from  [21,  Lemma  5]  we  have  that  the  first  ft  derivatives 
of  e  can  be  expressed  as 

e<,)  —  Ft  Vo,  CtH,  C,C,  e,)  ,  i=l,...,n,  (4.3) 


where  the  /i,’s  are  continuous  functions  of  their  arguments  and  the  e,’s  are  exponentially 
decaying  terms.  This  is  straightforward  to  show  (cf.  (2.48)-(2.49)  in  the  example  of  Sect.  2). 
starting  from  e  =  e  —  y0  and  using  (3.90)  and  the  facts  that  the  derivative  of  (  is  given  by 

n 

(  =  /4oC  +  we  -  P-1  J2  Cj MlClwwT Cj Mj Ci C  ,  ( 4.4 ) 

;=i 

and  that  the  elements  of  uy,  N,  and  £,  are  polynomial  functions  of  the  elements  of  C\H. 
Second,  from  (3.27)  we  have 


D‘ 

cJ  AqH  =  —  e>T(y,  u,  r)  +  e,(<) ,  0  <  i  <  n  -  1 , 

where  c,(<)  are  n^-dimensional  row  vectors  of  exponentially  decaying  terms.  Then. 
(3. 12)  to  express  o  as 


<PT(t/,u,r)  = 


Dt 1-2 


dth 


"{'.hr)  = 


On- 1  R  Dn_m_^D'  /x  Dn.n.l^Di  ,  , 

— r  -  — F - 2-j  P'°  (#)’•••»  — F - 2^*7=: — 


E  2 


El  1  ^  £22 


"21  ;=0^22 


being  bounded,  since  y  and  r  are  bounded.  Combining  (3.13),  (3.19)  and  (4.5)  we 
the  following  expressions  for  the  first  h  derivatives  of  the  tracking  error  e: 

bmD' 


(4.5) 

we  use 

(4.6) 

(4.7) 
obtain 


„(«)  _ 


Et 

bmD‘ 


[d>T(y, !/,  r )0  +  u]  +  e{t) 

p  [r/(y)u  +  d»T(y,u,r)(?]  +  e(<) 

bmD'  ,  bmEyD  D'  1  X,  \n  ,  ,,i 

— p—  [<7(.v M  +  — F - F  ^  (l/,  u.r)0  +  c(f) 

e-l 

bmD‘ 


E  i 


£i  Eo 
Eq  D  j  j  ,  _  [ 

£, 


[y(y)u]  +  :,-~CT,-i;,  lH  0  4-  e(f )  ,  1  <  >  <  n 


4.8) 
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It  is  important  to  note  that  bmEoD  / Ex  is  stable  and  proper.  The  boundedness  of  cT  A^H  is 
now  established  for  i  =  0, .  . . ,  fi  —  1  by  the  following  induction  argument: 

•  For  i  =  0,  (4.5)— (4.6)  give 

EjQ  h'Qll'2 

=  7T%'r)'  ~[q{y)u]  .  (4.9) 

£/o  ^0^2  ■U'l 

I  E  ~D 

Since  V(y,r),  —  (</(y)uj  are  bounded  and  — — £ —  is  a  row  of  stable  proper  filters.  (4.9) 
E\  EqE2 

implies  that  cT H  is  bounded.  Furthermore,  we  have  already  shown  that  e  is  bounded. 

•  For  1  <  i  <  n .  —  1,  assume  that  cT  AqH  and  are  bounded  for  0  <  k  <  i  —  1.  Hence. 
C,//  is  bounded,  and,  by  (4.3),  e*'*  is  bounded.  Then,  rewriting  (4.S)  as 


r  [<?(y)u]  =  T-e(,)  -  ^Tr~<?  A'0  lH0  + 


we  conclude  that  —  [</(t/)u]  is  bounded.  Finally,  using  (4.5)— (4  7),  we  obtain 
E  i 

ni  D'U 

cT  A'0H  =  -v(y,r),—-~[q{y)u} 

Esq  E'0^'2 

D' ^  E\Dn~2  D' 

=  —iy(y,r),———~[q(y)u]  +  t(t)  . 

E,  0  E'OE'2  Ej  i 

Hence,  cT  A'0H  is  bounded. 


(4.11) 


This  proves  that  H  is  bounded,  which,  by  (4.3),  means  that  e*'-,  1  <  i  <  ii,  are  bounded. 

Next,  we  prove  the  boundedness  of  u.  From  (3.13),  the  boundedness  of  0  and  the  fact  that 
U(y)  is  bounded  away  from  zero,  it  follows  that  u  is  bounded  if  u  and  4>(y,  u,r)  are  bounded. 
The  boundedness  of  II  implies  the  boundedness  of  p,  E  and  E_1 .  Since  y  =  =  (Cfi)-1E(’.',a, 

and  r.  and  C  are  bounded,  y  is  bounded  as  well.  Hence,  u  is  bounded.  From  (4.6) — ( 4 .7) . 
to  prove  boundedness  of  4>(y,u,r)  we  only  need  to  show  that  — pH -  [// ( ij)it}.  We  first  rewite 

E-2 

(4.8)  as 


-  lfl(y)«]  =  r~e<,)  ~  ~tt~cTTo  lH6  +  c(t),  1  <  i  s  n  -  l  - 


EoD  T  ,_1 


(4.1-2) 
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which  implies  that  —  [<?(y)u]  is  bounded  for  i  =  1, . . . ,  n  —  1.  Combining  this  with  the  fact 
Ci 

1 

that  it  follows  that  — -  [q(y)u]  is  bounded  for  i  =  l,...,n  —  3.  Differentiating 

Ci  C  2 

(4.12)  with  i  =  ft  —  1  and  using  (3.27)  and  (3.64)  we  obtain 


T-r  l?(yM  =  r-e(ri) 


^[cT/lofl//  +  ^(y,u,r)j0  +  e(0 


Substituting  or(y,u,  r)  from  (4.6)  and  rewriting  (3.11)  as  0T  =  [6»T  ,  ^On— 2]  1  we  express  (4.13) 


C  r  f  s  ,  DnEi  +  l 0n-2-Dn  ’Co  r  x  , 

rr[?(y)u)  =  - - mS/)“] 

^2  &1-C/2 

1  E  f)  ^ 

=  7  e<r1' - 7T-  cT^//0+  ^(y,r),  — —  [<?(yH  0  +e(C),  (4.14) 

^1  E2  J 

L  Dl 

which  implies  that  —  [<?(y)n]  is  bounded.  Since  L  is  of  degree  n  —  2  and  —  [</(y)n]  is 
C  2  C2 

Dn~ 2  jj  , 

bounded  for  i  =  1, ...  ,n  —  3,  it  follows  that  — — —  [q(y)u]  is  bounded.  Hence.  — C— 1  [r/(//)u] 

C2  Co 

is  bounded,  which  proves  that  u  is  bounded. 

In  order  to  show  the  boundedness  of  the  state  of  the  plant,  we  note  that  the  boundedness 
of  u  and  (3.16)  imply  that  D'y,  0  <  i  <  n  —  m,  are  bounded.  From  this  and  the  fact  that 
B(D)  is  Hurwitz,  we  conclude  that  the  state  x  in  (3.90)  is  bounded. 

We  have  thus  proved  that  the  state  of  the  closed-loop  adaptive  system  (3.90)  is  bounded 
on  [0,  <f).  Hence,  t(  =  00. 

To  prove  the  convergenceof  the  tracking  error  e  to  zero,  we  first  note  that  (3.39)  and  (3.47) 
imply  that  V  is  bounded  and  integrable  on  [0,  00).  Furthermore,  the  boundedness  of  t  (cf. 
(3.32)),  £  (cf.  (4.4))  and  H  (cf.  (3.90))  implies  that  V  is  bounded.  Hence,  V  — >  0  as  t  — ►  ec. 
which,  in  view  of  (3.47),  proves  that  e  — ►  0,  (  —< •  0  as  t  — >  00.  Since  77  =  5£  and  5  is 
bounded,  ij  —*  0  as  t  — ►  00.  Combined  with  (3.90)  and  the  boundedness  of  bm ,  this  also 
proves  that  r/0  — *•  0  as  t  —>  00.  Thus, 


lim  [y{t.)  -  yr(t. )]  lim  [e(t)  -  170(f)]  =  0 

t— *OG  t  —OO 


(4.15) 
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5  The  Class  of  Nonlinear  Systems 


Most  models  ol  nonlinear  systems  are  expressed  in  specific  state  coordinates.  From  that 
state-space  form  it  may  not  always  be  obvious  whether  or  not  the  nonlinear  system  at  hand 
has  the  input-output  description  assumed  in  Sect.  3.  Therefore,  we  now  give  coordinate- 
free  geometric  conditions  which  are  necessary  and  sufficient  for  a  single-input  single-output 
nonlinear  system  of  the  form 

z  =  f{z:a)  +  g(z\a)u  (ml) 

ij  =  h(z\a) 

to  have  an  input-output  description  of  the  form  (3.1),  which  is  repeated  here  for  convenience: 

m 

A{D)y  =  B(D)  [<?(t/)u]  +  Y,  D'  [p«o(y)  +  pj{y)0i\  •  (5-2) 

1  =  0 

In  (5.1 )  3  6  JRn  is  the  state,  u  G  1R  is  the  input,  y  £  R  is  the  output,  a  =  (oj  .  . .  arJ7  £  Hir 
is  a  vector  of  unknown  constant  parameters,  and  /,  g ,  h,  are  smooth  vector  fields  with 

/( 0;  a)  =  0,  h( 0;  a)  =  0,  for  all  a  £  Iff,  g(z)  ^  0  for  all  z  £  JR".  In  (5.2) 

•  the  coefficients  a0, . . .  ,an_l  of  the  denominator  polynomial  A(D )  =  Dn  +  an_lDn~y  + 

•  •  +  no  are  unknown, 

•  the  coefficients  60,  ...,bm(m  <  n  —  1)  of  the  numerator  polynomial  B(D)  ~  bmD’n  — 

•  +  b0  are  unknown, 

•  0 [  is  an  ^-dimensional  vector  of  unknown  parameters,  resulting  from  a  possible  overpa¬ 
rameterization  in  which  products  and  powers  of  the  original  unknown  parameters  n, 
are  treated  as  new  parameters  (so  that  i  >  r), 

•  P\j{y)'  0  <  i  <  m,  0  <  j  <  £  are  smooth  nonlinearities  with  q(y)  ^  0  Vy  6  IR. 

Pij( 0)  =  0,  0  <  i  <  m,  0  <  j  <  £. 

We  first  note  that  a  minimal  state  representation  of  (5.2)  is  given  by 

m 

x  =  Aex  +  bzq{y)u  +  X)rn-t  [p,-0(y)  +pJ(>j)Vi]  r 

1=0  l-t.4) 

T 

y  =  cl. v , 
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with  Av,  br,  cs,  as  defined  in  (3.91).  Hence,  the  following  statement  becomes  obvious: 

Fact  5.1.  The  nonlinear  system  (5.1)  has  an  input-output  description  of  the  form  (5.2)  if  and 
only  if  there  exists  a  global  in  z,  possibly  parameter-dependent,  diffeomorphism  transforming 
(5.1)  into  (5.3). 

Using  this  fact,  we  now  state  the  following  result: 

Proposition  5.2.  The  system  (5.1)  has  an  input-output  description  of  the  form  (5.2)  if 
and  only  if  the  following  conditions  are  satisfied  for  all  z  £  lRn  and  for  the  true  value  of  the 
parameter  vector  a: 


(Cl )  the  one-forms  dh ,  dl jh, . . .  ,  dLnj  lh  are  linearly  independent 


(C2)  [adyg,  adj</j  =  0  ,  t,j  =  0, . . .  ,n  —  1,  where  g  is  uniquely  defined  by 

0,  i  =  0 . n  —  2 


(C3) 


L9L)k  -  |  1?  i  =  n  _  1 


n  —  1 


adnfg  =  Y,  dj(a)ad}g  +  Y  [Pjo'iy)  +  V )  (y)0ij  ad/5 

j=0  j=o 

g,  ad Jfg\  =  0  ,  j  =  0, . . . ,  n  -  2 

m 

9  =  <?('/)HcJ(a)ad;/o^, 

j=u 

with  dj(o),Cj(a)  polynomial  functions  of  a.  &i  the  new  unknown  parameter  vector. 

fV 

a nd  pi(y)  —  /  p/l(v)dv.  i  =  0 . i 

Jo 

(Cl)  the  vector  fields  f  and  g  are  complete. 


Proof.  Using  Proposition  3  of  [23],  it  is  straightforward  to  show  that  conditions  (Cl)-((.3) 
are  necessary  and  sufficient  for  the  existence  of  a  local  diffeomorphism  such  that  in  the  new 
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coordinates  the  system  (5.1)  is  expressed  as 

(-l)rfB_i(a) 

I 

(-l)Vo(a)  0  ...  0 

0 

0 

(  —  1)"-™  [Prn0(y)  +  p£(j/)0i] 

(-i)n  [pooiy)  +  Po(y)^i] 

y  =  x  i , 

which  is  exactly  in  the  form  (5.3),  where  the  coefficients  a0, - a„_i,  b0, . . . ,  6™  depend  on 

the  physical  parameters  a.  From  [24],  condition  (C4)  is  necessary  and  sufficient  for  the  above 
diffeomorphism  to  be  global.  □ 

Remark  5.3.  The  above  proposition  gives  a  set  of  geometric  conditions  characterizing 
the  class  of  nonlinear  systems  to  which  our  adaptive  scheme  can  be  applied.  Whenever 
the  conditions  (0 i )-(C4)  can  be  verified  a  priori,  the  input-output  description  (5.21  of  the 
nonlinear  system  at  hand  is  determined  directly  from  (C3),  without  the  need  to  compute  the 
diffeomorphism  of  Proposition  5.2.  Unfortunately,  the  verification  of  (Cl )— ( C4 )  may  require 
some  a  priori  information  about  the  unknown  parameter  vector  a. 

Remark  5.4.  The  conditions  of  Proposition  5.2  are  satisfied  by  nonlinear  systems  that 
are  linearizable  by  output  injection  and  input-output  linearizable  by  full-state  feedback. 
However,  they  need  not  be  full-state  feedback  linearizable. 

We  illustrate  these  two  remarks  with  an  example. 
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Example  5.6.  It  may  not  be  obvious  that  (2.1)  is  the  input-output  description  of  the 
nonlinear  system 


~i  =  22  +  a  (3 yey  4-  2 y2) 

z2  =  23  -  Q  (2 yey  4-  y 2) 

23  =  24  4-  ayey 

~4  =  25  +  ay 2 

25  =  u  +  ay3 


(5.5) 


However,  this  can  be  established  by  checking  the  conditions  (Cl)-(C4).  Straightforward 
calculations  show  that  for  (5.5)  we  have 


y 

ad}  <7 
9 


•5-x - 4— - 1-  3— - 2— - h  -r — 

OZ\  OZ2  OZ3  OZ  4  OZ5 

a  [~3 y2g  +  (2 y  +  6y2)ad/y  -  {yey  +  ey  4-  4 y  +  3y2)ad}y] 
-g  4-  2ad fg  -  ad 2g  . 


(5.6) 

(5.7) 
(5.S) 


Hence,  the  conditions  (Cl)-(C4)  are  satisfied  for  all  a  and  the  input-output  description  of 
(5.5)  is 

D5y  =  (D2  +  2D  +  l)u  +  6[D2(yey  +  2y2  +  y3)  +  D(y2  +  2  y3)  +  y3j,  (5.9) 

where  0  =  a  It  is  important  to  note  that  to  determine  this  input-output  description  no 
explicit  change  of  coordinates  was  required.  In  this  simple  example,  however,  one  can  find 
the  corresponding  change  of  coordinates  by  inspection: 


-ci  =  2 1  +  2  z2  +  23 

X2  -  22  +  223  4-  24 

X3  -  23  4-  2 2.1  4-  25  ( •')■  10 ) 

x.j  =  24  4-  22.5 

X  5  —  . 
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In  these  coordinates.  (5.5)  becomes 


X\  —  X’2 

±2  =  2  3 

£3  =  £4  +  u  +  a  ^ yey  -j-  2 y1  -f-  y3)  (5.11) 

£4  =  £5  +  Q  ( y 2  +  2y3) 

£5  =  u  +  ay 3 

y  =  xi. 

We  immediately  see  that  (5.11)  has  the  input-output  description  (5  9). 

However,  it  should  also  be  pointed  out  that  (5.5)  is  not  full-state  feedback  linearizable. 
since  the  distribution 

Q3  =  span  jg,  ad/g,  ad^g,  ad3g  j  (5.12) 

is  not  involutive[3,4j.  □ 

6  Conclusions 

This  paper  has  extended  the  theory  of  adaptive  control  for  linear  systems  to  a  class  of  systems 
which  are  essentially  nonlinear  in  the  sense  that  their  nonlinearities  are  not  restricted  by  any 
growth  constraints.  In  spite  of  this  absence  of  growth  constraints,  all  the  stability  and 
tracking  results  are  global. 

The  assumptions  on  the  linear  part  of  the  system  are  the  same  as  in  the  standard  adaptive 
theory  for  linear  systems.  However,  to  guarantee  the  aforementioned  global  properties,  the 
systematic  design  procedure  has  departed  from  the  two  main  ingredients  of  most  adaptive 
schemes  for  linear  systems:  the  certainty-equivalence  control  and  the  normalization  of  the 
update  law.  In  addition  to  the  certainty-equivalence  part,  the  control  contains  a  term  which 
counteracts  the  effects  of  rapidly  growing  nonlinearities.  Thanks  to  the  presence  of  this 
term,  the  normalization  of  the  update  law  is  avoided,  which  allows  the  rapid  decrease  of 
the  parameter  error.  This  proved  to  be  crucial  in  preventing  finite  escape  times  common  in 
systems  with  rapidly  growing  nonlinearities. 


The  class  of  nonlinear  systems  has  been  restricted  bv  coordinate Tree  geometric  conditions 
which  are  equivalent  to  the  structural  requirements  that  the  nonlinearities  depend  only  on 
the  output  and  do  not  enter  the  system  before  the  control  does,  and  that  the  zero  dynamics 
are  linear  and  exponentially  stable.  Relaxing  these  restrictions,  and  thus  enlarging  the  class 
of  systems  that  can  be  adaptively  controlled  using  only  output  measurement,  is  a  topic  oi 
further  research. 
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